More Evidence for a Distribution of Tunnel Splittings in Mni 2 -acetate 
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In magnetic fields applied parallel to the anisotropy axis, the magnetization of Mni2 has been 
measured in response to a field that is swept back and forth across the resonances corresponding 
to steps N — 4, 5, ...9. The fraction of molecules remaining in the metastable well after each sweep 
through the resonance is inconsistent with expectations for an ensemble of identical molecules. The 
data are consistent instead with the presence of a broad distribution of tunnel splittings. A very 
good fit is obtained for a Gaussian distribution of the second-order anisotropy tunneling parameter 
Xe = — In (| E | /2D). We show that dipolar shuffling is a negligible effect which cannot explain 
our data. 
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PACS numbers: 75.50.Xx, 75.45. +j 

Much attention has been focussed on the behavior of 
high-spin molecular nanomagnets because of their in- 
trinsic interest as well as for their potential applica- 
tions in high density storage of information and pos- 
sible use as qubits for quantum computation. In 
the prototypical single molecule magnet Mni2-acetate, 
([Mn 12 Oi 2 (CH 3 COO)i 6 (H 2 0)4]-2CH3COOH-4H 2 0), re- 
cent studies [1] of the magnetic relaxation at temper- 
atures below 0.5 K for different sweep rates of a mag- 
netic field applied along the longitudinal (easy) axis have 
yielded evidence that the symmetry breaking field that 
drives the tunneling in this material is locally varying, 
broadly distributed second-order anisotropy which is for- 
bidden by the tetragonal symmetry of the Mn 12 -acetate 
crystal, but which is present in real crystals due per- 
haps to long-range crystal dislocations or local chemical 
environment. The distributed second order transverse 
anisotropy gives rise to a distribution of tunnel splittings. 
In this paper, we present additional experimental evi- 
dence that the tunnel splittings are broadly distributed. 

Mni 2 -acetate is composed of weakly-interacting mag- 
netic clusters, each consisting of twelve Mn atoms cou- 
pled by superexchange through oxygen bridges to give 
a net spin S = 10 at low temperatures. The magnetic 
molecules are regularly arranged on a tetragonal body- 
centered lattice with strong uniaxial anisotropy of the or- 
der of 60 K. Below the blocking temperature, Tb ~ 3 K, 
steep steps are observed [2] in the M versus H curves due 
to enhanced relaxation of the magnetization whenever 
levels on opposite sides of the anisotropy barrier coincide 
in energy. The steps are separated by plateaux where 



magnetic relaxation requires thermal activation and is 
strongly suppressed. 

The Hamiltonian for Mn 12 -acetate is given by: 



H 



-DS 2 Z - ASt - g z n B H z S z + V T , 



(1) 



where the anisotropy D = 0.548(3)lf, the fourth-order 
longitudinal anisotropy A = 1.173(4) x 10~ 3 K, g z is es- 
timated to be 1.94(1) [3], and Vt contains all symmetry- 
breaking terms responsible for the observed tunneling. 
For temperatures below approximately 0.5 K, it has been 
shown that essentially all the tunneling proceeds from the 
ground state level m' = —10 of the metastable well [4]. 
Ground state tunneling, indicated by the arrow in the 
doublewell structure in Fig. 1, occurs at magnetic fields: 
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where N = — (m + m') is the step number. The data 
reported in this paper were obtained at 0.25 K, where 
only ground state tunneling is observed. 

The magnetization of small single crystals of Mni 2 - 
acetate was determined by methods described elsewhere 
[5]. The relaxation of the magnetization was studied for 
each observable ground state resonance as follows: start- 
ing from zero magnetization, the external field was swept 
from zero up to and beyond a given resonance N. When 
the field reached the middle of the plateau just above 
step N, it was reversed and swept back across this reso- 
nance. When the field reached the middle of the previous 
plateau, the field was swept back up and the process was 
repeated many times. The same procedure was repeated 
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for all observed resonances, N = 4,5, ...9, as shown in the 
normalized magnetization curves of Fig. 1. The turning 
points were chosen so that the amplitude of the field oscil- 
lations were the same for each TV and the turning points 
of the field occurred near the center of a plateau. The 
field was swept at a rate dH z /dt = 5.83mT/s. 



consistent with the behavior expected for a sample com- 
prised of an ensemble of identical molecules. 
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FIG. 1. Normalized magnetization versus magnetic 
field for a field that is swept back and forth across 
each resonance. The field was swept at a rate of 
dH z /dt = 5.83mT/s. The upper inset shows the 
labeling described in the text for the first few os- 
cillations for N = 7. The arrow in the lower inset 
designates tunneling from the ground state of the 
metastable well for step N = 5. 

As the field is swept back and forth across the field res- 
onance, the magnetization relaxes in steps. The manner 
in which it relaxes is expected to be different for a set of 
identical molecules than for the case of a distribution of 
tunnel splittings. We will show below that an ensemble 
of identical molecules cannot describe the data presented 
here. A method of fitting different distributions to the 
magnetization data will be developed, and an excellent fit 
is obtained for a lognormal Gaussian distribution. It will 
also be argued that dipolar shuffling [6] would produce 
an asymmetry in the steps which is not observed. 

The fraction of molecules remaining after each step can 
be determined from the magnetization curves. At the 
turning points labeled by fo, fx, fi,—, for N — 7, shown 
in the upper inset of Fig. 1, the fraction of molecules 
remaining in the metastable well is related to the nor- 
malized magnetization by, 



(3) 



where Mj is the normalized magnetization at the turning 
points of the field. For each N, Fig. 2 shows the fraction 
of molecules remaining in the metastable well after each 
field pass plotted as a function of how many times, j, the 
field was swept past the resonance. The issue is whether 
or not the fraction remaining after each field sweep is 
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FIG. 2. Fraction of molecules remaining in the 
metastable well as a function of the oscillation num- 
ber, j. With a Gaussian distribution of Xe = — ln(| 
E | /2D), the fits to Eq. 5 for each N shown by 
the solid lines are nearly indistinguishable from the 
measured data. The upper inset shows the first few 
data points for N — 6 and 7 on a linear scale to 
illustrate that dipolar shuffling has minimal effects 
(see text). The lower inset shows the dependence of 
the fitting parameter, Xq, on 7Y. 

For a set of identical molecules, the probability, P/v, of 
remaining in the metastable well for resonance, N, is the 
same for each molecule. If the fraction of molecules in 
the metastable well before sweeping past the resonance 
is /at.o, then the fraction of molecules remaining in the 
metastable well after sweeping the field through a reso- 
nance will be /jy,i = fN,oPN- When the field is swept 
back down across the resonance a second time, another 
set of molecules will tunnel. The fraction remaining after 
the second pass can be determined from the fraction re- 
maining after the first pass, fx,2 = fN,iPN = /at,o(-P/v) 2 - 
In general, the fraction of molecules remaining in the 
metastable well after j passes of the field is: 



fN,j = fN,o( p N) J 



(4) 



This form implies that for a set of identical molecules, 
fft,j would depend exponentially on j, and a semi-log 
plot of the remaining fraction, /atj, versus the number 
of oscillations, j, should yield straight lines for each N. 
As can be seen in Fig. 3, which shows the data of Fig. 2 
on a semi-logarithmic scale, this is not the case. The 
simple assumption that all the molecules are identical is 
clearly inconsistent with our data. 
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FIG. 3. The fraction of molecules remaining in 
the metastable well should have an exponential de- 
pendence on the oscillation number, j, for an en- 
semble of identical molecules (see text). Since the 
data is not linear on a semilogarithmic plot, our re- 
sults are inconsistent with expectations for identical 
molecules. 

We now consider the case of a distribution of tunnel 
splittings. Chudnovsky and Garanin [7] calculated the ef- 
fect of such a distribution due to transverse anisotropics 
caused by crystal dislocations which produce deforma- 
tions that require inclusion of a magnetoelastic term in 
the Hamiltonian. The dislocations give rise to a trans- 
verse anisotropy, E(S% — Sy) that varies with the distance 
from the dislocation. Thus, the second order transverse 
anisotropy prefactor, E, is distributed. Alternatively, 
Cornia et al. [8] have pointed out that there are six differ- 
ent isomers of Mni2-acetate, four of which are only two- 
fold symmetric. This lowers the symmetry group for the 
Hamiltonian and also permits inclusion of locally varying 
second-order transverse anisotropy terms. In either case, 
a distribution of second-order anisotropy gives rise to a 
distribution of tunnel splittings. The following analysis 
is valid regardless of the origin of such a distribution of 
the second order transverse anisotropy prefactor, E. 

The probability of remaining in the metastable well 
is given by the incoherent [9] Landau-Zener-Stueckelberg 
formula [10], P/v = exp(— irAff /2hv]y), where A/v is the 
tunnel splitting and vn = (25 — N)g/iBdH z /dt is the en- 
ergy sweep rate. The tunnel splitting can be calculated 
through high order degenerate perturbation theory. For 
the case of second order transverse anisotropy (neglect- 
ing the fourth order longitudinal anisotropy), the tun- 



nel splitting has the form, A/v = ?1n9n 



where 



(instead of D = 0.548 K) to compensate for the neglect 
of the fourth order term. For even N, £ N = S — N/2 and 
r) N = 1. For odd N,£ N =S-(N- l)/2 and rj N = N/2. 

For an arbitrary initial distribution, fo(E), the frac- 
tion of molecules remaining in the metastable well with 
a transverse anisotropy prefactor between E and E + dE 
is h{E)dE = P 1 (E)f (E)dE, where P X (E) is the 
probability of remaining in the metastable well for the 
N = 1 resonance. Integrating over all possible trans- 
verse anisotropy prefactors, one obtains the fraction of 
molecules remaining in the metastable well after hav- 
ing swept the field entirely through the resonance, fi — 
J-™ Pi(E)fo{E) dE. The total fraction of molecules re- 
maining in the metastable well after having swept the 
field back and forth across the N = 1 resonance j times 
is, fi,j = (Pi(E)Y fo(E) dE. Note that for N > 1, 
the field is swept once through each of the lower reso- 
nances. Therefore, the fraction of molecules remaining 
in the metastable well, after sweeping the field back and 
forth j times across resonance iV is given by, 

h,j = / II Pk(E) (Pn(E)Y f (E) dE, (5) 
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where rifclo 1 Ps(-^-e) determines how much of the distri- 
bution remains after having swept through the previous 
N— 1 resonances and P/v (E) is the probability of remain- 
ing in the metastable well for the N th resonance. Eq. 5 
can be used to model the data shown in Fig. 2. Different 
distributions, fo(E), can be chosen and all the curves in 
Fig. 2 should be fit simultaneously for the same choice of 
parameters defining the shape of fo (E) . 

An excellent fit to the data is obtained for a log-normal 
distribution of the form [1] f(X E ) = e -(x E -x„) 2 /4a 2 ^ 
where Xe = — ln( E \ /2D). This form captures the be- 
havior remarkably well for all observed resonances, and 
is indistinguishable from the data shown in Figs. 2 and 3. 
However, rather than being constant, the parameter Xq 
required for a best fit was found to vary weakly with step 
number N. The values of Xq chosen for each TV in the fit 
are shown in the inset to Fig. 2. The width of the distri- 
bution, (7 = 0.09814, was constrained to be constant for 
the fit. The variation of Xq with TV could be due perhaps 
to the existence of other symmetry-breaking terms in the 
Hamiltonian. In our previous studies [1], deviations from 
full scaling were attributed to a small admixture of tun- 
neling due to transverse magnetic fields. 

"Dipolar shuffling" has been recently proposed as an 
explanation for an unexpected dependence of tunnel 
splittings on sweep rate [6]. This effect refers to the fact 
that when a molecule tunnels, the internal dipole fields 
are locally altered so that some molecules that were pre- 
viously below are now above the resonance field and vice- 
versa. For the long thin sample studied here (a w 40/im, 
b w 40/im, c « 150/im), this is expected to introduce an 



9n = [(2S-n-2)\^ \] {2S (N)\ 2SV ' and we set D = °- 65 K 
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asymmetry between up- and down-sweeps: during an up- 
sweep many of the dipole-shuffled molecules will skip the 
resonant field resulting in a lower tunneling rate while 
during a down-sweep molecules will have additional op- 
portunities to tunnel giving an enhanced tunneling rate. 
This would yield a tilted staircase behavior in the frac- 
tion of molecules with successive values of j. The upper 
inset of Fig. 2 shows that the fraction remaining varies 
smoothly with oscillation number and does not exhibit 
a tilted staircase behavior. We conclude that for the 
samples studied here, there was no observable effect as- 
sociated with dipolar shuffling, which cannot therefore 
be responsible for the nonexponential behavior shown in 
Fig 3. 

In summary, the dependence of the measured tun- 
neling probability on the amount of magnetization re- 
maining in the metastable well was studied by sweeping 
the field back and forth across each field resonance N. 
The fraction of molecules remaining in the metastable 
well after each sweep through the resonance was shown 
to be inconsistent with expectations for an ensemble of 
identical molecules. The data were analyzed on the as- 
sumption that there is a distribution of tunnel splittings 
associated with locally varying second-order transverse 
anisotropy arising from dislocations or a local symmetry- 
breaking term due to detailed chemical environment. A 
very good fit is obtained for a Gaussian distribution in 
X E = - In (| E | 1 2D). A variation of X a on N is found 
which may be due to the existence of other symmetry- 
breaking terms in the Hamiltonian. In principle, the 
analysis presented here offers a method for determining 
the distribution of tunnel splittings. 
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